   
   cap program drop SFA2tier    
   SFA2tier ln二手房 ln财政收入 ln收入 ln赤字  ln城市投资 ln销售面积 ln住宅投资 ln医院 ln人口密度 lngdp ln绿化, search     
   version 10.0
   ml max, difficult 
   est store m1     //模型3
   
      SFA2tier ln二手房 ln财政收入 ln收入 ln赤字  ln城市投资 ln销售面积 ln住宅投资 ln医院 ln人口密度 lngdp ln绿化  ///
            序号1 序号2 序号3 序号4 序号5 序号6 序号7 序号8 序号9 序号10 序号11 序号12 序号13 序号14 序号15 ///
			序号16 序号17 序号18 序号19 序号20 序号21 序号22 序号23 序号24 序号25 序号26 序号27 ///
			序号28  序号29  序号30 ///
            , search   
   version 10.0
   ml max, difficult 
   est store m2    //模型4
   
   
      SFA2tier ln二手房 ln财政收入 ln收入 ln赤字 ln城市投资  ln销售面积 ln住宅投资 ln医院   ///
            ln人口密度 lngdp  ln绿化  年份2007-年份2021            ///
            , search   
   version 10.0
   ml max, difficult 
   est store m3  //模型5
   
   
      SFA2tier ln二手房 ln财政收入 ln收入 ln赤字 ln城市投资 ln销售面积 ln住宅投资 ln医院   ///
            ln人口密度 lngdp ln绿化   年份2007-年份2021            ///
			序号1  序号2  序号3 序号4 序号5 序号6  序号7 序号8  序号9 序号12 序号13 序号14  序号16 序号17  ///
			 序号20   序号25  序号27     ///
            , search   
   version 10.0
   ml max, difficult 
   est store m4  
  
   SFA2tier_sigs 


 *---------------  
*-估计单边效应 
*---------------  
  dropvars *_hat* *_diff* pat doc doc_pat
  SFA2tier_eff
  replace uw_diff_exp = -uw_diff_exp 
  
  foreach v of varlist u_hat w_hat uw_diff *_exp{
     replace `v' = `v'*100
  }
 
   *-样本总体
     tabstat u_hat_exp w_hat_exp uw_diff_exp,     ///
             s(mean sd p25 p50 p75) f(%6.2f) c(s) save
    tabstatmat all
    mat all = all'
    mat list all, title(full sample)
	
    
   *-年度统计
      bysort 年份: tabstat u_hat_exp w_hat_exp uw_diff_exp, ///
             s(mean sd p25 p50 p75) f(%6.2f) c(s) 
       tabstat  uw_diff_exp, by(年份) ///
             s(mean sd p25 p50 p75) f(%6.2f) c(s)    
  
                        
   *-urban: 东中西
     bysort 地区: tabstat u_hat_exp w_hat_exp uw_diff_exp, /// 
             s(mean sd p25 p50 p75) f(%6.2f) c(s)   

			 
			 
 *-直方图  u的图
   set scheme s1mono
   histogram u_hat_exp, percent title(Percent, place(10) size(*0.7)) ///
       ylabel(,angle(0)) ytitle("") xtitle("供给者 (%)")  ///
       xscale(titlegap(3) outergap(-2))   
       graph export Result\u_his.wmf, replace

   
   histogram w_hat_exp, percent title(Percent, place(10) size(*0.7)) ///
       ylabel(,angle(0)) ytitle("") xtitle("需求者 (%)")  ///
       xscale(titlegap(3) outergap(-2)) 
       graph export Result\w_his.wmf, replace   /// w的图
	   
	   
	   * replace uw_diff_exp = -uw_diff_exp
   histogram uw_diff_exp, percent title(Percent, place(10) size(*0.7)) ///
       ylabel(,angle(0)) ytitle("") xtitle("净剩余 (%)")   ///
       xscale(titlegap(3) outergap(-2))    
       graph export Result\uw_his.wmf, replace  /// 图3
	   
	   